Correlated few-photon transport in one-dimensional waveguides: linear and nonlinear 
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We address correlated few-photon transport in one-dimensional waveguides coupled to a two- 
level system (TLS), such as an atom or a quantum dot. We derive exactly the single-photon 
and two-photon current (transmission) for linear and nonlinear (tight-binding sinusoidal) energy- 
momentum dispersion relations of photons in the waveguides and compare the results for the different 
dispersions. A large enhancement of the two-photon current for the sinusoidal dispersion has been 
seen at a certain transition energy of the TLS away from the single-photon resonances. 

PACS numbers: : 42.79.Gn, 42.50.-p, 03.65.Nk, 32.80.Qk 



I. INTRODUCTION 

The study of equilibrium dynamics of correlated pho- 
tons and mixed photonic-atomic excitations (called po- 
laritons) has received huge interest in the last few years. 
It allows one to simulate strongly correlated condensed 
matter phenomena, for example, quantum phase tran- 
sition in simple controllable quantum systems [1], such 
as coupled optical cavities with each containing a single 
two-level system (TLS) [2]. It is especially interesting to 
model nonequilibrium dynamics in these systems. Now, 
it is a highly nontrivial task to study out-of-equilibrium 
dynamics of many correlated photons nonperturbatively. 
Instead, we pose a relatively simpler problem: we exam- 
ine the exact dynamics of few correlated photons in one 
dimension. 

The dynamics of few correlated photons in the Dicke- 
type Hamiltonian with a tight-binding (TB) sinusoidal 
energy dispersion relation for photons was investigated 
recently using a computational technique [3] . The single- 
photon dynamics in that model was studied before in 
Ref. [4] to demonstrate a quantum switch for the coherent 
transport of a single photon. The two-photon scattering 
states in a similar model for a linear energy dispersion 
relation of photons was derived exactly by employing 
a generalized Bethe ansatz method [5] and Lehmann- 
Symanzik-Zimmermann reduction techniques [12]. Re- 
cently large-scale ultrahigh-Q coupled resonator arrays 
with a TB sinusoidal dispersion relation have been real- 
ized in photonic crystals [7]. 

Here we study correlated few-photon transport for a 
TB sinusoidal dispersion relation with band edges using 
an analytical method [8, 9]. We then compare the re- 
sults for the nonlinear sinusoidal dispersion relation with 
a linear dispersion relation. We apply the Bethe ansatz 
method to study the linear dispersion [5, 10, 11]. There 
are some studies on correlated photon transport for co- 
sine dispersion [12], but none has tried before to compare 
correlated dynamics in the waveguides for different dis- 
persions. Shi, Fan, and Sun [13] considered correlated 
two-photon transport in a linear waveguide coupled to 
cavity with a TLS to study a photon blockade in the 
strong TLS-cavity coupling regime [14]. It is possible to 



further extend the method of this paper to investigate 
correlated photon dynamics in a cavity coupled to a non- 
linear waveguide. 



II. MODEL 

We consider a system consisting of two one- 
dimensional (ID) coupled-resonator-optical waveguides 
being connected by a TLS, such as an atom or a quan- 
tum dot [4]. The two-photon or multiphoton dynamics 
in this system is strongly correlated [3, 5, 11]. The full 
Hamiltonian within the random-phase approximation is 
given by 

oo 

Tin = -J ^2 '(a\.a x+ i_ + al +1 a x ) +W.a z /2 
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Here, a\ (a x ) denotes the photon creation (annihilation) 
operator at site x, and J is the strength of hopping be- 
tween nearest-neighbor sites. ^Z' indicates summation 
over all integers, omitting x = — 1, 0. The Pauli operators 
a z and a± (= a x ± ia y ) represent the TLS at the zeroth 
site. The TLS with a transition frequency fl is coupled to 
the left and right waveguides by Vl and Vr, respectively. 
Here we assume a direct coupling between the waveguides 
and the TLS. A generalization of the present results for 
a sidccoupling is straightforward. The energy dispersion 
of photons in the waveguides is given by Ek = — 2 J cos k 
with a wave number — n < k < n. Henceforth, subscript 
or superscript n and / denote properties pertaining to 
nonlinear and linear dispersions, respectively. 

Now, we redefine in terms of a new Hamiltonian rl, 
where we replace the TLS by an appropriate additional 
bosonic system at the zeroth site [3, 15]. Then rl = 
Ho + V, with 

OO 

Ho = -J (dla x +i + a] B+1 a at ) + W$b 
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+ Vb(oLi6 + 6 t o_i) + V a ((46 + 6 t oi) , 
V = | 6+6(6^-1) , (2.2) 
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where b is the photon annihilation operator at the ze- 
roth site. The ground and excited states of the TLS in 
H n correspond to no photons and a single photon, re- 
spectively, at the zeroth site of H. The unwanted mul- 
tiphoton occupancies at the TLS have been avoided by 
introducing the interaction term V. We expect that the 
Hamiltonian in Eq.(2.2) in the limit U — > oo is equivalent 
to the Hamiltonian H n . We define the current operator, 
i x = —i V x (a\.a x +i — a x+1 a x ), using a continuity equa- 
tion, where V x = — Vl (— Vr) for x = — 1 (1), and V x — J 
for all other x. 

A linear energy dispersion for photons, i.e., E k — v g k 
is a good approximation in many situations for ID opti- 
cal waveguides, such as microwave transmission lines [16] 
and surface plasmon modes of a metallic nanowires [17]. 
Here, v g is the momentum-independent group velocity of 
photons. The real-space Hamiltonian of the ID waveg- 
uides and TLS for a linear energy dispersion is given by 



Hi 



-IV Q 



E 

CB=1,2 



(x) + htttfb 



-(V lC {(0)b + V 2 cl(0)b- 



u, 



H.c.) + -tfb(tfb- 1) (2.3) 

where c\(x) [c\(x)] is a bosonic operator creating a pho- 
ton at x in the left-hand [right-hand] side of the impu- 
rity. We define the current operator in the system as 
1 = -«[%, iVi - N 2 ]/2, where N t [N 2 ] is the total num- 
ber of photons in the left-hand [right-hand] side of the 
TLS. Then we derive I = i[Vic\(0)b- F 2 c|(0)6- H.c.}/2. 
One can readily check that both current operators I x 
and I are equivalent. We need to derive exact scatter- 
ing eigenstates of Hi to find an average steady-state cur- 
rent in the system. Employing a standard transformation 
to even-odd modes, c\{x) — \V\c e (x) + V 2 c (x)]/V and 
c 2 (x) = [V 2 c e (x) - ViCo(x)]/V with V = VW+Vi, the 
Hamiltonian in Eq.(2.3) breaks into two decoupled parts, 
i.e., Hi — H e + H , where 



H e 



U 
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dx c\{x)—c e {x) + Mlb% 

tfb(tfb-l) 
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F(4(0)6 + 6t Ce ( )), 

d_ 

dx' 



and 



dx c\{x)jr-c {x). 



III. SINGLE-PHOTON DYNAMICS 



It is simple to find the single-photon scattering state 
for the Hamiltonian in Eq.(2.2). Let us denote the single- 
photon scattering state by 4>k(x) (= (x\(f)k}) at lattice 
site x with incident momentum k. For a photon being 
incident from the left (a right-moving photon), i.e., < 



k < 7T, we find 4>k{x) — e lkx 
qe lkx for x > 0, and fe (O) 



f" 

T-k- 



2iV L V R smk 



-r%e lkx for x < 0, ^fc(z) 
■■-J{l+rl)/V L with 



e ik (vZ + v£)-j(nn-E k y 



v, 



■t n k - 1 (3.1) 
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FIG. 1: (Color online) Plot of the single-photon current 
(I n (k),I l (k)) with Q for a linear (L) and a tight-binding (TB) 
energy dispersion with incident photon energy (a) = 1.01 
and (b) Ek = 1.6. Here, Et, fi, V , and V" are in units of J, 
with J — 1 and h = 1. 



The transmission coefficient, TJ} = |i£| 2 shows a Brcit- 



Wigner-like (i.e., Lorentzian) line shape around the res- 
onance transition energy HQ — Ek(2J 2 — V% — V R l )/2J 2 . 
For symmetric coupling, i.e., Vl = Vr = V , a single 
photon at the single-photon resonance can fully trans- 
mit from one side of the impurity to the other side. The 
single-photon current for the nonlinear sinusoidal disper- 
sion is given by I n (k) = (4>k\I x \4>k) = 2JT£ sinfc, where 
2Jsinfc is the velocity of a photon. We plot I n (k) for 
different V and Ek in Fig.l. 

The single-photon scattering state |1, fc) of Hi for an 
incoming photon from the left with momentum k is given 

by 



dx 



Vi 

, } {y[9k(x)cl(x) 



V 



6(x)e k tf] + -±h k (x)cUx)}\0,0), 



V 



where g k (x) = e lkx [9{-x) + T k 9{x)], h k {x) = e lkx ,e k = 
V/(vgk-hQ + iV 2 /2v g ), r k = e k /e* k , with 6(x) being the 
step function. \n, m) denotes n photons in the waveguide 
and to photons at the impurity site (TLS). We determine 
the transmission coefficient T l k for a photon from the left 
to the right lead by rearranging the single-photon state 
in terms of the original field operators. It is given by 



rpl 
1 h 



4r!r 2 



{v g k-hn) 2 + (r 1 + r 2 ) 2 



(3.2) 



(2.4) with r a = V 2 /2v g . The single-photon reflection coeffi- 
cient is R{ = (V^ 4 + V 2 4 )/V 4 + 2V 1 2 V 2 2 Re[T fc ]/V 4 . 

We can also compute T l k from Eq.(3.1) after linearizing 
the energy dispersion along with a proper regularization 
scheme. This occurs at the matching condition for the 
wavevector Ik ~ 7r/2; thus, we find E k ~ ±2Jfc — Jtt 



by expanding the cosine around k = ±7r/2 and v g = 2J 
(and lattice constant I = 1 here). We further require 
renormalization of the coupling as Vl = Vi/2 and Vr — 
V 2 /2 in Eq.(3.1). The steady-state single-photon current 
for the linear dispersion is given by I (k) = (1, fc|/|l, k) — 
VgTf./2-K. We need to multiply I l (k) by 2n to compare 
I (k) with I n (k). This is due to our convention for the 
normalization of the scattering states. The current I l (k) 
also shows a Breit-Wigner-type line shape around the 
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FIG. 2: (Color online) Plot of the real (ReKo) and the imag- 
inary (lmKo) parts of _Ks k (0) for E kl = 1.01, i?fc 2 = 1.6, and 
V' — 0.2, 0.3. Inset shows a magnification of the rightmost 
part of the plot. E k2 , ^, and V' axe in units of J, with 

J = 1 and h = 1. 
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FIG. 3: (Color online) Plot of the real (ReTfo) and the imag- 
inary (Imifo) parts of KE k (0) for two sets of parameters. 
Sl= {E kl = 1.01, E k2 = 1.6}, S2= {E kl = 1.6, E k2 = 1.6} 
with V' = 0.2 in both sets. E kl , E k2 , fi, and V' are in units 
of J, with J — 1 and ft = 1. 



resonance Ml — v g k. We plot I l {k) (after multiplying by 
2tt) for different V x = V 2 = V" and E k in Fig. 1 and 
compare it with that for the TB sinusoidal dispersion. 
We find that the line-shapes of the single-photon current 
are qualitatively similar for both linear and TB sinusoidal 
dispersions. The differences in the width and the height 
of the current line shapes for the two dispersions are due 
to the variance in the group velocity, which is strongly 
fc-dependent for the sinusoidal dispersion. 



IV. TWO-PHOTON DYNAMICS 

Tight-binding sinusoidal dispersion: We study two- 
photon transport in the lattice model using a tech- 
nique based on the Lippman-Schwinger scattering theory. 
This technique was developed recently to discuss electron 
transport in quantum impurity models [8, 9]. We here 
assume that the zeroth site is initially empty. The two- 
photon scattering states 'J'k(x) of the Hamiltonian H for 
arbitrary value of U are given by (see Appendix A) 

where x = (xi,x 2 ), k = (ki,k 2 ), = (0,0), the 
symmetrized initial state </>k(x) = [4>k 1 {xi)(j)k 2 {x2) + 
4>ki ( x 2) ( f>k 2 {xi)]/V2, and the total energy of the two in- 
cident photons is E^ = — 2 J(cos k\ +cos k 2 ). The expres- 
sion A^E k (x) = (x|Gj(i?k)|0) has the explicit form 



w d qi dg 2 q (x)0*(O) 
2j_ n J_ n (2tt)2 E^-E^+ir, 



We find the two-photon scattering state ^(x) of the 
Hamiltonian H n by taking the limit U — > oo in Eq. (4.1): 

M*) = (x|^> = k (x) - ■ (4.2) 



We can immediately see from Eq.(4.2) that V'k(O) = 0; 
i.e., both photons cannot be at the TLS (or zeroth site) 
simultaneously, but the occupation probability of a sin- 
gle photon at the zeroth site has a finite value. A single 
photon at the zeroth site corresponds to the excited state 
of the TLS. The average occupation of the impurity site 
is calculated by taking the expectation value of the op- 
erator n b (= 6+6) in |^ k ). (n b ) = £ X1 \{x u 0|Vk)| 2 . The 
nature of the scattering states V'k(x) depends on KE k (0). 
Interestingly, we find that both the real and the imagi- 
nary parts of KE k (0) become zero for a broad range of 
the transition energy Ml away from any of the single- 
photon resonances; for example, KE k (0) jumps to zero 
near Q, = 1.92 for V = 0.2 as shown in Fig. 2. We plot 
the real and the imaginary part of KE k (0) in Figs. 2 and 
3 for different incident energy and coupling strength. The 
zero of KE k (0) depends on the strength of coupling V' 
(see Fig. 2) but is independent of the incident energy, 
as shown in Fig. 3. We also notice that the real part of 
KE k (0) passes through zero between two single-photon 
resonance energies. 

The two-photon current (or transmission) has been cal- 
culated exactly by taking the expectation value of the 
current operator I x in the scattering state \ipu.) = |0k) + 
|5k), where |Sk) = Gg (Ek)V\ipk) is the interaction- 
induced correction to the scattering state. The two- 
photon current I n (ki,k 2 ) (= (V'kl^lV'k)) has two parts: 
one, j™(ki,k 2 ) (= (^klXrl^k)), is the contribution from 
two noninteracting photons, and the other, 5j n {k\,k 2 ) (= 
3 c + js)i is induced by the nonlinear photon-photon 
interaction at the impurity site. Thus, I n (ki,k 2 ) = 

+ 5j n (ki,k 2 ). Here, Jc (= (0 k |4|S k ) + 
(Skl-fzl^k)) i s an expectation value of the current op- 
erator between scattered and incident photon wave func- 
tions, and it is a measure of two-photon cross correlation. 
An expectation value of I x in the interaction-induced 
scattered wave function is given by js = (Skl^rlSk)- We 
now determine all these terms separately. We find that 
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FIG. 4: (Color online) Plot of the two-photon current (transmission) change due to interaction at the TLS site for a linear (L) 
and a tight-binding (TB) sinusoidal dispersion. The energies of the incident photons are as follows: (a) E kl ~ 1.01, E k2 = 1.6, 
(b) E kl — 1.01, E k2 — 1.6, and (c) E kl — 1.6, E k2 — 1.6. E kl , E k2 , O,, V' , and V" are in units of J, with J = 1 and h=l. 



jf(k u k 2 ) = Af[I n (ki) + I n (k 2 )}, where TV is the total 
number of sites in the system. For i > 1 or 2 < -1, we 
find j s = 2 J Im (5 k |at a x+ i|Sk) with 

m I t lq \ 2|^ k (o)| 2 r dg [M 2 

(5 k |4a x+ i|S' k ) = — / —1^(0)1 I (q)I 1 (q), 



where / m (t?) 



* dqi ^(0)fc(^™) 



2tt £^-£; n 



17] 



, 772 = 0,1, 



and, j c = 2 J Im (0 k |(aj.a x 



+i 



L x+1 L 



k |4 i a :E2 |S' k ) 



^ 2 (^l)C g (°) 



IjSk) with 



^(^1)^,(0) 



77/ 



Eh, — E„ 



77/ 



Thus, jf(ki, k 2 ) is one order of magnitude higher than 
Sj n (k±, k 2 ) for a finite system [8]. In the redefined system 
described by the Hamiltonian in Eq. (2.2), the two pho- 
tons see each other only when both are at the impurity 
site, i.e., at x — 0, and the amplitude for this to occur 
is on the order of 1/Af for a finite system of AT sites. In 
Fig. 4, we plot the interaction-induced two-photon cur- 
rent change 6j n (ki,k 2 ) for different values of V and 
energy of the incident photons. The current change 
5j n (ki, k 2 ) is nonzero near the single-photon resonances. 
We find that the direct correlation js has the same sign 
(which is positive for transmitted photons) for all values 
of il, whereas the cross correlation jc can have different 
signs depending on f2. A two-photon bound state results 
in the two-photon scattering state for both the linear [5] 
and the nonlinear dispersions [8, 18]. The bound state 
acts as a composite object and remains together after 
passing through the TLS. It has a maximum contribution 
in the two-particle current change for Ek t = Ek 2 = Ml, 
which is manifested by js in Fig. 5. The energy and mo- 
mentum of the scattered photons are redistributed over a 
wide range of values satisfying the total energy conserva- 
tion [19]. This redistribution of energy and momentum 



emerges as the background fluorescence [5], which can 
be conceived as being a result of the inelastic scattering 
of one photon from a composite transient object formed 
by the TLS and the other photon [20]. The magnitude 
of 5j n (ki, k 2 ) around the single-photon resonance is quite 
large when the energy of the incident photons is the same 
[see Fig. 4(c)]. We find in Fig. 5 that the cross correlation 
jc has a large negative value near the single-photon res- 
onance. A photon spends more time near the resonance 
and thus the effective strength of photon-photon interac- 
tion is increased much when both the incident photons 
have the same energy. One can also interpret the large 
reduction of 5j n (ki, k 2 ) near the single-photon resonance 
as being due to antibunching of the scattered photons [5] 
or as a manifestation of the photon blockade. 

Surprisingly, we find a large enhancement of 
&j n (ki,k 2 ) for a certain value of the transition energy M7, 
which is away from any of the single-photon resonances. 
It occurs at the value of Ml where KE k (0) becomes zero, 
and it is relatively more pronounced for the incident pho- 
tons with different energies as shown in Figs. 4(a) and 
4(b). This is a special "two-photon resonance" arising for 
a finite bandwidth of the dispersion [9, 21]. The direct 
correlation term js only contributes in 8j n [k\, k 2 ) at this 
special two-photon resonance. We can also phrase that a 
special "two-photon bound state" is formed for the finite 
bandwidth TB dispersion, and it creates a strong bunch- 
ing of the transmitted photons. Interestingly, the magni- 
tude of Sj n (ki,k 2 ) at the two-photon resonance depends 
on the position of the nearest single-photon resonance. 

Linear dispersion: Now we determine the two-photon 
dynamics of Hi- It was recently studied [5, 11] using 
the Bethe ansatz approach for a linearized dispersion of 
the Hamiltonian H n - We here construct the two-photon 
scattering eigenstates for Hi following Ref. [11] for an 
arbitrary U, and we show that our results in the limit 
U — » oo indeed match with the earlier results of Refs. 
[5, 11]. The general two-photon scattering eigenstate is 
of the form: 
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dx\dx 2 



A 2 {g{x l7 x 2 )-^=cl(x 1 )cl(x 2 ) + e(x 1 )S(x 2 )cl(x 1 )b'' + w S(x 1 )5(x 2 )-^=b^b' i } + B 2 {t(x 1 ;x 2 )c< e (x 1 )cl(x 2 ) 



1 



+f(x 1 )5(x 2 )cl(x 1 )tf}+C 2 h(x 1 ,x 2 )^cl(x 1 )cUx 2 )\\0,0) 



(4.3) 



r 




Sj (k,, k 2 ) 
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FIG. 5: (Color online) Two-photon current change j$, jc, an d 
5j n (ki, fe) at x > for a tight-binding sinusoidal dispersion 
with energy of the incident photons Ek ± ~ Ek 2 = 1.6 and 
V' = 0.2. Ekx, Ek 2 , O and V' are in units of J, with J = 1 
and ft = 1. 



where g(x\,x 2 ) and h(xi,x 2 ) are the probability ampli- 
tudes for both the photons in the even or the odd sub- 
space, respectively, while t(x\,x 2 ) is the amplitude for 
one photon in the even and the other in the odd sub- 
space. The quantity e(x) (f(x)) is the probability am- 
plitude for one photon in the e (o) subspace when the 
impurity is occupied by another photon; w is the prob- 
ability amplitude for both the photons at the impurity 
site. Here, A 2 , B 2 and C 2 keep track of the incident pho- 
tons. We determine the amplitudes in Eq. (4.3) using the 
two-photon Schrodinger equation. The results are given 
in Appendix B. We find, as U — > oo, w — > 0, and the 
other amplitudes of the scattering state in Eq. (4.3) re- 
duces to that of the TLS- waveguide obtained in Ref. [11] 
with the original Dicke-type Hamiltonian. Thus, we find 
an indirect proof for the validity of the mapping scheme 
for the TLS-waveguides on the Anderson impurity model 
with infinite U in the present nonequilibrium dynamics 
study. 

Here we estimate the two-photon current by taking 
the expectation value of I between the state in Eq. 
(4.3) in the limit U oo. We choose both the in- 
coming photons from the left of the TLS; i.e., A 2 — 
V 1 2 /V 2 ,B 2 = V X V 2 /V 2 , and C 2 = V 2 2 /V 2 . We find 
I l (k u k 2 ) = C[l\ki) + /%)]/2tt + 6j l (h,k 2 ), where C 
is the length of the full system, and 



Sj l (h,k 2 ) 



YiYi 

ir 2 V 5 
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2 _* 



(4.4) 



to photon-photon interaction created by the TLS for the 
linear dispersion. We plot Sj l (ki,k 2 ) (after multiplying 
by 4-7T 2 ) in Fig. 4 with the transition energy Ml for dif- 
ferent coupling strengths and incident energies. We find 
that 5j l (ki,k 2 ) is qualitatively similar to 8j n {k\,k 2 ) near 
the single-photon resonances. We again notice a large re- 
duction of the two-photon current near the single-photon 
resonance when incident energy of the two photons is 
the same [see Fig. 4(c)]. The current change 5j l (Ek 1 — 
E k2 = W) = -l6V?V 2 2 vl/(n 2 V s ) is independent of the 
incoming photon's energy. The physical mechanism of 
the large reduction of Sj l (ki,k 2 ) near the single-photon 
resonance is similar to that of the TB dispersions. How- 
ever, Sj l (ki,k 2 ) is almost zero near Ml ~ 2 J for the 
linear dispersion. Thus there is no special two-photon 
resonance away from the single-photon resonances for the 
linear dispersion without band edges. It is quite differ- 
ent from the two-photon current line shape for the TB 
sinusoidal dispersion. The quantitative distinction be- 
tween 6j l (ki,k 2 ) and 6j n (ki,k 2 ) near the single-photon 
resonances can be attributed again to the variance in the 
group velocity of a photon for the linear and TB sinu- 
soidal dispersions. 



where 5j l (ki,k 2 ) is the two-photon current change due 



CONCLUDING REMARKS 



In summary, we have shown qualitative and quanti- 
tative differences in the correlated two-photon transport 
in ID waveguides for linear and TB sinusoidal disper- 
sions. A special "two-photon resonance" in the two- 
photon transmission line shape of the sinusoidal disper- 
sion has been found. It might be useful to achieve a 
strong photon-photon interaction for nonmonochromatic 
light at low intensity. There have been enormous ad- 
vances in experiments to realize optical nonlinearities at 
the few-photon level by creating strong light-matter in- 
teractions in various systems [16, 17, 22]. We hope that 
our prediction of two-photon resonance for a TB sinu- 
soidal dispersion will be observed in the near future. We 
also plan to study the correlated dynamics of few pho- 
tons in a cavity quantum electrodynamics setup with a 
two-level or multilevel (driven) atom (or atoms) and pho- 
tons with tight-binding sinusoidal dispersion, such as, in 
one-dimensional coupled-resonator-optical waveguides. 
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Appendix A: Derivation of the state in Eq. (4.1) 

The full scattering eigenstate of the Hamiltonian W = 
Ho + V is given by the Lippman-Schwinger equation 

\1>) = \<l>)+G$(E k )V\1>), (Al) 

where the noninteracting two-photon Green's function is 
Go(Eis.) = l/(Ek — Ho +iv)- I R the position basis |x) (= 
a xi a x 2 1 )' wnere I ) is a vacuum), we obtain 

Vk(x) = k (x) + UK Ek (x) MO), (A2) 

where Vk(0) = </>k(0)/[l - UK Ek (0)} is derived from Eq. 
(A2) by inserting x = 0. 



Appendix B: Amplitudes of the state in Eq. (4.3) 

The amplitudes in Eq.4.3 are as follows: 

9(xi,x 2 ) = ^=9^1)9^X2) - -Jj-CV £ ^^ 

x e - l (hn- l v*/2v g )( Xl - X2 )/v ge ( Xi _ X2 )Q( X2 ) + (1 ^ 2), 

e{x) = ^{gk 1 {x)e k2 + g k2 {x)e kl ^j 

+ C' e t( - Ek ~ hn+ty2 / 2v ^ x / Vg 9(x) 
t(xr,x 2 ) = {g kl {xi)h k2 (x 2 ) + g k2 (x 1 )h kl (x 2 ))/2TT, 

and f(x) = (e kl h k2 (x) + e k2 h kl (x))/2ir, h(x 1 ,x 2 ) = 

(h kl (xi)h k2 (x 2 ) + h k2 (xi)h kl (x 2 ))/2ir\/2, where C = 

(iV/TTV g )e kl e k2 + f3(e kl + e k2 )/(n(l - /?)) with (3 = 
{iV 2 /2v g )/{W. + U/2 - Ek/2). 

t Present address: Department of Physics, University 
of Cincinnati, Cincinnati, Ohio 45221, USA 
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